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ABSTRACT 

We use new and archival Chandra and ROSAT data to study the time variability of the X-ray 
emission from the pulsar wind nebula (PWN) powered by PSR B1509—58 on timescales of one week 
to twelve years. There is variability in the size, number, and brightness of compact knots appearing 
within 20" of the pulsar, with at least one knot showing a possible outflow velocity of ~ 0.6c (assuming 
a distance to the source of 5.2 kpc). The transient nature of these knots may indicate that they are 
produced by turbulence in the flows surrounding the pulsar. A previously identified prominent jet 
extending 12 pc to the southeast of the pulsar increased in brightness by 30% over 9 years; apparent 
outflow of material along this jet is observed with a velocity of ~ 0.5c. However, outflow alone cannot 
account for the changes in the jet on such short timescales. Magnetohydrodynamic sausage or kink 
instabilities are feasible explanations for the jet variability with timescale of ~ 1.3—2 years. An arc 
structure, located 30" — 45" north of the pulsar, shows transverse structural variations and appears to 
have moved inward with a velocity of ~ 0.03c over three years. The overall structure and brightness 
of the diffuse PWN exterior to this arc and excluding the jet has remained the same over the twelve 
year span. The photon indices of the diffuse PWN and possibly the jet steepen with increasing radius, 
likely indicating synchrotron cooling at X-ray energies. 

Subject headings: ISM: individual (G320.4—1.2) — ISM: jets and outflows — pulsars: individual 
(B1509—58) — stars: neutron — supernova remnants — X-rays: ISM 


1. INTRODUCTION 

The pulsar B1509—58 and the supernova remnant 
(SNR) G320.4—1.2 (MSH 15—5^) represent one of ap¬ 
proximately 20 known associations between a pulsar and 
a SNR. This pulsar is one of the most energetic known, 
with a period (P) of 150 ms, a period derivative (P) 
of 1.2x10“^^ s s“^, a characteristic age Tc = P/2P w 
1700 yr, a spin-down luminosity E = 47r^JPP“^ w 1.8 x 
1037 g,.gg g-i moment of inertia / = 10"^® g cm^), 

and an inferred dipole surface magnetic field Bp w 
3.2 X lO^^iPpy/^ » 1.5 X 1013 G llKasni et allll^ 
iLivingstone et al]l2005ll . 

SNR G320.4—1.2 has been well studied at radio, op¬ 
tical, and X-ray wavelengths. The radio morphol¬ 
ogy consists of a partial shell to the southeast and 
a series of bright c lumps ~ 25' to the northwest 
l|Gaens1er et af.l Il999t) that coincide with the optical 
nebula RCW 89 liRodgers. Campbell, fc WhiteoaUflOBnl: 

iSeward et alJll983il. The X-ray morphology consists of a 

bright, elongated pulsar wind nebula (PW N) with a colli¬ 
mate d jet extending ~ 4' t o the s o utheast llSewai det_alJ 
1983t (Greiveldin^ Tet^ 119951: [T russoni et all Il9^ 
Brazier fc Becked Il997t IGaensler et al.1 12002L hereafter 

^ Alfred P. Sloan Research Fellow 


G02). To the northwest are therm al clumps associate d 
with the rad io clu mps and RCW 89 llSeward et alJl98^ . 
In addition, IG02I identified several compact knots close 
to the pulsar, plus two semicircular arcs at a distance 
of 17" and 30" to the north of the pulsar. The toroidal 
morphology and collimated jet are reminiscent of struc¬ 
tures found i n the PWNe pow e red by the Crab and 
Vela pulsars llHe ster et all 119951 i Weisskopf et aLlEoonl: 
iHelfa.nd. Gotthelf. V Ha,lr)ernll2flfll|l . The arcs have been 
proposed as being due to ion-compression in the particle- 
dominated equatorial flow from the pulsar and were in¬ 


terpreted bv lG02l as analogs of the “wisps” found in the 
Crab Nebula. SNR G320.4—1.2 has recently been ob¬ 
served in very high energy y-rays by HESS. The emis¬ 
sion is elongated along the PWN axis possibly indicat- 
ing inverse Compton s cattering of relativistic electrons 
4,ronia,n et a,lJ 


The structures in the Crab and Vela PWNe are known 
to vary in brightness and position over short timescales 
(days to months). In the case of the Crab Nebula, the 
outward moving X-ray and optical wisps (with velocity 
V ~ 0.5c) are thought to mark the PWN termination 
shock, while the small-scale X-ray and optical knots are 
thought to ident ify unstable, quasi -stationary shocks in 
the pulsar wind llHester et al.ll200'^ . Radio wisps, which 
rarely correspond to optical wisps, develop and move out- 
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ward at slower velocities {v ~ 0.3c) and there are even 
more slowly moving radio f eatures (v ^ 10'^ km s~ ^) far¬ 
ther away from the pulsar llBietenholz et alJl2004ll . Re¬ 
cent observations of the wisps and polar knots in the near 
infrared indicate bri ghtness variations on time scales as 
short as 20 minutes ijMelatos et al.ll2?)f)^ . In the case of 
the Vela PWN, the X-ray arcs also move outward and 
vary in brightness by up to 30% llPavlov et al.ll200lll . 

The variability of the Vela PWN jet observed in X- 
rays has been attributed to both kink instabilities, to 
account for the dramatic shape and brightness changes 
over the course of days, and sausage instabilities, to ac- 
count for the relat ivistically moving “blobs” {v ^ 0.5c; 
iPavlov et al.ll200^ . The Crab nebula jet, on the other 
hand, shows only weak X- ray morphologica l variations 
on year-long time scales lIMori et all 1200411 and rela- 
tivistic outflow ide ntified in the optical with v ~ 0.4c 
llHester et al.ll20d^ . 

The PWN powered by PSR B1509—58 represents a 
unique opportunity to study variability in PWNe. The 
physical size of the PWN is approximately 10 and 100 
times larger than the Crab and Vela PWNe, respectively. 
Given the observed variability in the Crab and Vela 
PWNe and the distance to PSR B1509—58 (5.2±1.4 kpc; 
iGaensler et al.ll99^ . lG02l predicted measurable variabil¬ 
ity in the PWN of PSR B1509—58 on time scales of a 
few years. To that end, we obtained new Chandra X- 
ray Observatory observations of PSR B1509—58 and its 
PWN. In this paper we compare our new images to exist¬ 
ing ROSAT PSPC and HRI images and Chandra ACIS-I 
images and report on the variations observed over time 
scales of I week to 12 years. 

2. X-RAY OBSERVATIONS AND ANALYSIS 

We observed G320.4—1.2 with the Chandra AGIS-I de¬ 
tector on 2003 Apr 21 for a total of 9.6 ks, on 2003 Apr 28 
for a total of 10 ks, and on 2003 Oct 18 for a total of 
19.4 ks. The data were calibrated using CIAO version 3.1 
and GALDB version 2.28. After filtering for good time 
intervals, the total useable exposure times were 9469 s, 
9497 s, and 18739 s for the Apr 21, Apr 28, and Oct 18 ob¬ 
servations, respectively. Exposure corrected images were 
constru cted f or each epoch between 0.3 and 8 keV as out¬ 
lined in lG02L The telescope was not dithered during the 
2003 Apr 28 observation so that the effects of the chip 
gaps and dead columns could not be removed at that 
epoch. 

G320.4—1.2 had been previously observed with the 
Chandr a AC IS-I detector on 2000 Aug 14 for a total of 
19.3 ks IIG02I1 . The archival data products were obtained 
from the Chandra archive and were recalibrated using 
CIAO version 3.1 and GALDB version 2.28 to make use 
of the newest gain solutions and correct the geometry of 
the data set^. The data were filtered for good time inter¬ 
vals resulting in a final exposure time of 17863 s and an 
exposure corrected image was made as outlined above. 

We compare these Chandra data to archi val obser¬ 
vatio n s with the Rdntqen Sat ellite (R OS A T: iTriimneil 
Il982t iPfeffermann et all Il987r . G320.4—1.2 was ob¬ 

served with the ROSAT PSPC detector on 1991 Feb 22 

^ For standard ACIS data preparation, see 

http: / / cxc. harvard.edu / ciao / guides / acis_data. html 
For intormation on the geometry correction, see 

http: / / cxc.harvard.edu/cal/Hrma/optaxis / platescale 


- Mar 8 and 1992 Feb 25 - 27 for a total o f 9.1 ks 
jGreiveldinger et al.jll9^ iTrussoni et~m 1199611. Data 
were also taken with the ROSAT HRI detector on 
1994 Feb 10 - 23 and 1^ 94 Sep 15 for a total of 22.5 ks 
ijBrazier fc BeckeJ Il997jl . The archival data products 
were obtained from NASA’s HEASARC^. Exposure cor¬ 
rected images were made using the Xselect utility and 
the hriexpmap and pcexpmap tasks within version 5.3.1 of 
the FTOOLS package. The exposure times of the PSPC 
and HRI data are approximately matched to produce the 
same number of counts on the two images. To facilitate 
comparison to the Chandra images, exposure corrected 
images were also made from 2 ks each of the 2000 Aug 14 
and 2003 Oct 18 Chandra data sets to match the energy 
range (~ 0.1 — 2.5 keV) and number of counts^ on the 
ROSAT PSPC and HRI images. 

The four Chandra images were registered to each other 
using a brightness-weighted mean of PSR B1509—58, the 
star Muzzio 10 (located at (J20 00) R.A. 15^13™55!2, 
deck -59°07'51"6: IMiizzinlll^ . and four other back¬ 
ground sources in the field (located at (J2000) R.A. 
15'^13™41®, deck -59°11'45"; R.A. 15^14^00®, deck 
-59°12'38"; R.A. 15^14’^05®, deck -59°14'40"; R.A. 
15'^14™32®, deck -59°08'09"). The relative regis¬ 
tration error is < 07l with rotation constrained to 
< 1°. Although the pulsar does not have a significant 
proper motion (un < 39 mas yr“^, yis < 52 mas yr“^, 
IGaensler et aklflO^ . its suitability for registration pur¬ 
poses is in question because the strong pile-up enhances 
the azimuthal brightness asymmetry that results from a 
misali gnment of the tel e scope mirrors in th e innermost 
shell (iJarius et akll200nl: IPavlov et akll2001ll . Therefore, 
we also determined the registration solution excluding 
the pulsar. Again using a brightness-weighted mean, we 
find a relative registration error < 071 for the 2000 Aug 
14 and 2003 Oct 18 images and < 0725 for the two 2003 
Apr images. 

Spectra were extracted for selected sources of inter¬ 
est using the FTOOLS Xselect utility for the ROSAT 
PSPC data and the CIAO acisspec script for the Chan¬ 
dra data. Background spectra were extracted from an 
annular region surrounding each source of interest. The 
spectra were rebinned to have a minimum of 20 counts 
per spectral channel. XSPEC version 11.3.1 was used for 
the spectral analysis. Unless otherwise noted, we present 
the 90% confidence limit for errors. 

All physical size and dist ance calculations assu me that 
the source is 5.2 kpc away l|Caensler et aklll9^ . 

3. SUMMED CHANDRA IMAGE 

Shown in Figure ^ is a summed and exposure cor¬ 
rected image of 50 ks of the Chandra data excluding the 
2003 Apr 28 observation, which was not dithered. The 
energy range is 0.3—8.0 keV and the image has been con¬ 
volved with a Gaussian of FWHM 10". The pulsar po¬ 
sition is indicated by the white cross. In addition to the 
jet and inner and outer arcs (G02 features “C”, “5”, and 

http://heasarc.gsfc.nasa.gov/ 

Version 3.6a of the PIMMS program 
|http://cxc.harvard.edu/toolkit/pimms.jsp I was used to de¬ 
termine the Chandra exposure time needed to match the number 
of counts on the ROSAT image for the diffuse PWN and knots in 
the RCW 89 region. 








































3 


-1-'-r 


-59 00 - 





RCW 89 


03 - 


o 

o 

o 

CM 


U 

u 

Q 


06 - 


arc b 


arc a 


\ 


09' - outer arc 


12 - 



inner arc 


•59 15 1=_ L 


15 " 14 "^ 40 ^ 


20 ® 


X 

14 ^ 10 ^ 


40® 


20 ® 


X 

13“0® 


RA (J2000) 


Fig. 1.— Summed 50 ks Chandra image of G320.4—1.2 over the energy range 0.3—8.0 keV. The image has been exposure-corrected and 
convolved with a Gaussian of FWHM 10^^ The transfer function is linear. On this and subsequent images, the pulsar position is indicated 
by the white cross. Overlaid on the jet are the regions used to determine the radial variation of the photon index as indicated in Table ^ 


“E”, respectively), we identify two faint arc-like struc¬ 
tures at distances of ~ 110" and 160" from the pulsar. 
These newly identified features are labelled as “arc a” 
and “arc b” and are ^20" in width. Arc a can be iden¬ 
tified in the individual 20 ks 2000 Aug and 2003 Oct 
images and also appears in more recent Chandra images 
(P. Slane, private communication). The newer Chandra 
data will require more analysis before we can determine 
if arc b is also a robust structure. 

4. TIME EVOLUTION 
4.1. Flux of the Diffuse PWN 

To determine if the total flux of the diffuse PWN 
varied over time, we performed simultaneous absorbed 
power law fits to the ROSAT PSPC and the Chandra 
data sets, requiring the same absorbing column {Nh) 
and photon index (P) for all data sets, but allowing the 
normalization to vary between epochs (to account for 
both calibration uncertainties and actual, evolutionary 


brightness variations). We define the diffuse PWN to 
be the low brightness emission surrounding the pulsar, 
interior to the RCW 89 region to the northwest, and 
extending slightly beyond the jet to the southeast. We 
specifically exclude the jet and inner and outer arc struc¬ 
tures from this definition. The exact region u sed fo r the 
analysis is the same as shown in Figure 7 of IC02L The 
best fit values are Nh = (8.6 ± 0.2) x 10^^ cm“^ and 
P = 1.97 ± 0.03. Tabled shows the unabsorbed fluxes 
obtained from the simultaneous spectral fit for the en¬ 
ergy range 0.5—10.0 keV.^ Note that the loss of dither¬ 
ing for the 2003 Apr 28 observation resulted in lower 
flux values for the diffuse PWN and jet. The total unab¬ 
sorbed flux of the diffuse PWN has remained steady at 
« 5.8 X 10“^^ erg cm“^ s“^ from 1991/1992 to 2003. 

® Although the ROSAT PSPC data do not extend above 2.5 keV, 
fluxes were extrapolated up to 10 keV using a model spectrum with 
the same photon index and normalization found from the spectral 
fits. 
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TABLE 1 

Spectral Fits to Various Subregions of the Source 


Region 

s/p 

R/C‘> 

Epoch 

r 

Fx^ 

(10“12 erg cm-2 g-l) 

Diffuse PWN 

s 

R 

1991/1992 

1.97±0.03 

57±3 


s 

C 

2000 Aug 14 


58±3 


s 

C 

2003 Apr 21 


58±3 


s 

c 

2003 Apr 28=^ 


53±3 


s 

c 

2003 Oct 18 


58±3 


I 

c 

2000 Aug 14 

1.97±0.03 

58±3 


I 

c 

2000 Apr 21 

1.98±0.03 

59±3 


I 

c 

2003 Apr 28"^ 

1.96±0.03 

53±3 


I 

c 

2003 Oct 18 

1.96±0.03 

58±3 

Jet 

s 

R 

1991/1992 

1.64±0.07 

4.1±1.0 


s 

C 

2000 Aug 14 


5.5±0.6 


s 

C 

2003 Apr 21 


5.2±1.0 


s 

c 

2003 Apr 28'i 


4.1±1.0 


s 

c 

2003 Oct 18 


5.4±0.6 


I 

c 

2000 Aug 14 

1.62±0.11 

5.5±0.6 


I 

c 

2003 Apr 21 

1.53±0.11 

5.3±1.0 


I 

c 

2003 Apr 28'i 

1.79±0.13 

4.0±1.0 


I 

c 

2003 Oct 18 

1.65±0.11 

5.4±0.6 

Outer Arc 

s 

R 

1991/1992 

1.64±0.09 

4.4±1.0 


s 

C 

2000 Aug 14 


3.7±0.6 


s 

C 

2003 Apr 21 


3.4±1.0 


s 

C 

2003 Apr 28 


3.3±1.0 


s 

c 

2003 Oct 18 


3.5±0.6 


Note. — Uncertainties are all at 90% confidence. The models used are power 
law of the form cc E~^ where P is the photon index and the integrated flux 
is Fx = F{E1, E2) = Ef{E)dE. All models assume interstellar absorption 
using the cross sections of lBalucihska-Church &: McCammonI 119921) . assuming solar 
abundances. 

^“S” indicates the results of simultaneous fits to multi-epoch data requiring the 
same P, holding Nh fixed at 8. 6 x 1 0^^ cm~^ and allowing the normalization to vary 
between epochs as outlined in ini “I” indicates a fit to an individual epoch holding 
Njj fixed at 8.6 x 10^^ cm“^. 

'^“R” indicates ROSAT observations, “C” indicates Chandra observations. 

^Fluxes are for the energy range 0.5—10.0 keV, and have been corrected for inter¬ 
stellar absorption. 

“^The telescope was not dithered during the 2003 Apr 28 observation resulting in 
a decrease in measured flux for those structures that were intersected by chip gaps 
and dead columns. 


The value of the absorbing column does not match 
that of 10021 within the errors. We fitt ed just the recal¬ 
ibrated 2000 data set studied by IG02I and found Njj = 
(8.2 ± 0.3) X 10^^ cni“^ and T = 1.91 ± 0.04, which is 
consistent with the multi-epoch fit. We also f ixed Nr 
at 9.5 X 10^^ cm“^ and T at 2.05 (as found bv 10021) to 
determine what, if any difference there is in total flux. 
For the energy range 0.5—10.0 keV, we obtained an un¬ 
absorbed flux of 5.8 X 10“^^ erg cm“^ s“^ for both sets 
of values for the absorbing column and the photon index 
indicating that these parameters are sufficiently degener¬ 
ate that the slightly different values do not significantly 
affect the fi tted fl ux. This flux is within the error of that 
reported bv 10021 We conclude that the mismatch in our 
abs orbin g column and photon index values compared to 
the 10021 results is most likely due to updated Chandra 
calibration solutions. Unless otherwise indicated, spec¬ 
tral fits in this paper hold Nh fixed at the simultaneous 
epoch fit value of 8.6x 10^^ cm“^. Fluxes are reported for 
an energy range of 0.5—10.0 keV and have been corrected 
for interstellar absorption. 

4.2. Flux of the Pulsar 


To determine if the pulsar flux has varied over time, 
we compared the background subtracted count rate from 
a 1'.'4 radius region around the pulsar for each Chan¬ 
dra epoch over the energy ranges 0.5—10.0 keV and 
2.0—10.0 keV. Note that these count rates do not re¬ 
flect the true flux from the pulsar because of pile-up. We 
chose the two different energy ranges to determine if any 
variation in count rate could be due to pile-up differences 
between epochs. A reduction in pile-up might result from 
the extra absorption at low energies caused by the time- 
dependent ACIS filter contamination. The count rates 
are shown in Table |21 In 2000, the pulsar count rates 
between 0.5—10.0 keV and between 2.0—10.0 keV were 
~8% less than in 2003, representing less than a ia differ¬ 
ence. Because the difference is approximately the same 
for the two energy ranges, it is unlikely that we have a 
significant bias from a reduction of pile-up at low ener¬ 
gies. We also determined the flux by performing a power 
law fit to each of the Chandra dat a sets correc ting for in¬ 
terstellar absorption and pile-up l|Davyi200llf . We held 
Nh fixed at 8.6 x 10^^ ^ fixed at 1.19 

llOusumano et al 1 IM1 for each epoch but allowed the 
normalization and pile-up grade morphing parameter to 
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vary. As shown in Table [21 these fits indicate that the 
flux of the pulsar has remained steady at ~ 5.1 x 10“^^ 
erg cm“^ s“^. 


4.3. Small-Scale Structure Near the Pulsar 
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Fig. 2.— Chandra images showing changes in the environment 
immediately surrounding the pulsar. The 2003 Apr image combines 
both April observations. The images have been convolved with a 
Gaussian of FWHM 1". The energy range is 0.3—8.0 keV. The 
grayscale range is linear and ranges from 0 to 3% of the peak value 
on all three images. The black arc indicates the location of the 
outer edge of the inner arc (see Figure [Tl. Arrows on 2003 images 
indicate a possible knot in motion. 


Figure [2I shows changes in the structure immediately 
surrounding the pulsar between 2000 Aug and 2003 Oct®. 
In 2000, the structure consisted of four compact knots < 
3" across (corresponding to a physical size of 0.1—0.2 pc) 
and located along the jet axis between distances of 
~ 3" — 17" on both sides of the pulsar. Three years 
later in 2003 Apr, this structure has completely changed, 
with approximately nine small, unresolved knots (< 0.5", 
< O.Olpc) again located along the jet axis between dis¬ 
tances of ~ 3" — 17" on both sides of the pulsar. No 
apparent changes were evident between 2003 Apr 21 
and 2003 Apr 28. Approximately six months later in 
2003 Oct, the structure has again changed, with six knots 
ranging in size from 0'.'5 to 4" and again located along the 


® Note that the small-scale structure identified here was beyond 
the resolution capability of the ROSAT detectors. 


jet axis within 17" of the pulsar. In all cases, the knot 
activity is greater to the northwest than to the southeast 
and in no case are knots found beyond the inner arc at 
a radius of 17". 

There is the possibility that some of the knot variabil¬ 
ity we observe could be due to artifacts. It is known that 
artificial structures may appear near sources with strong 
pile-up. However, we believe that most of the variabil¬ 
ity is real because the typical count rates are ~ 0.005- 
0.05 counts s“^ (corresponding to 100-1000 counts on 
the 2000 Aug and 2003 Oct images) - higher than would 
be expected for artificial structures. Another indication 
that the small-scale knots are genuine comes from our 
Monte Carlo simulations discussed in (ESI In each of 
our ten Poisson noise simulations, prior to smoothing 
to 10", we recover all of the knot structures. Further¬ 
more, in analyzing their knot 1, G02 specifically ruled 
out asymmetries in the wings of the point-spread func¬ 
tion and pulsar photons assigned to the wrong location 
on the sky as possible explanations. Finally, the small- 
scale knots are apparent on unsmoothed images and are 
not due to the 1" FWHM Gaussian convolution applied 
to Figure 121 

The knot structure changes so drastically that identify¬ 
ing motion is practically impossible. However, if the knot 
to the southeast of the pulsar in 2003 Apr and 2003 Oct 
is indeed the same feature (identified with arrows on Fig¬ 
ure O, then the 4" motion results in a velocity of ~ 0.6c 
(assuming outflow along the jet with a 30° inclination^ 
to the line of sight and correcting for relativistic Doppler 
boosting). 


4.4. The Jet 

We performed simultaneous absorbed power law fits to 
the ROSAT PSPC and Chandra data fo r the jet in the 
same manner as for the diffuse PWN in ii4.ll We fixed 
Nh at 8.6 X 10^^ cm“^ and required the same photon 
index for all the data sets but allowed the normalization 
to vary between epochs. The photon i ndice s and fluxes 
are indicated in Tabled As noted bv lG02L the photon 
index of the jet is flatter than the diffuse PWN. From the 
spectral fits, we note that the jet brightened by ~30% be¬ 
tween the ROSAT PSPG observation in 1991/1992 and 
the Chandra observations, however the flux remained 
steady from 2000 to 2003. The flux difference between 
the 1991/1992 ROSAT data and the 2000 Chandra data 
is slightly greater than the 90% confidence limit. This 
flux difference might be caused by the differing spectral 
responses of ROSAT and Chandra, however this is pre¬ 
sumably accounted for by the corresponding auxiliary 
response files. We can also exclude a systematic uncer¬ 
tainty in the calibration of the two observatories as con¬ 
tributing to the flux difference by examining the flux of 
the Diffuse PWN. As seen in Tabled the flux from the 
extended emission remains constant as expected. 

We also compared the 1991/1992 ROSAT PSPG 
image, 1994 ROSAT HRI image and 2000 Aug and 
2003 Oct Chandra images of the jet as shown in Fig¬ 
ure 121 To facilitate comparison between different ob- 

^ |G03 derive a 30° jet inclination to the line-of-sight based 
on Doppler boost ing and radio polarization arguments, while 
lYatsu et alJ I20Q5I’) derive an inclination angle > 50° based on an 
interaction between the unseen northwest jet and the RCW 89 re¬ 
gion. 
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TABLE 2 

Detected Count Rate and Flux for PSR B1509-58 


Epoch 

Total Count Rate^ 
(counts s~^) 

Hard Count Rate^ 
(counts s“^) 

Flux'= 

(10“^^ erg cm~^ 

2000 Aug 14 

0.130±0.005 

0.103±0.005 

5.0±1.8 

2003 Apr 21 

0.144±0.008 

0.115±0.008 

5.1±1.8 

2003 Apr 28 

0.138±0.008 

0.108±0.008 

5.1±1.8 

2003 Oct 18 

0.142±0.005 

0.111±0.005 

5.2±1.8 


^Count rates are for the energy range 0.5—10.0 keV. They have been background 
subtracted but have not been corrected for pile-up. Errors are from Poisson statis¬ 
tics. 

^Count rates are for the energy range 2.0—10.0 keV. They have been background 
subtracted but have not been corrected for pile-up. Errors are from Poisson statis¬ 
tics. 

^Fluxes are for the energy range 0.5—10.0 keV, and have been corrected for in¬ 
terstellar absorption and pile-up. Uncertainties are all at 90% confidence. The 
models used are power law of t he form /r oc E ~^ wh ere P is the photon index 
and has been held fixed at 1.19 iCusumano et aT 112001 and the integrated flux is 
Fx = F(E1, E2) = Ef{E)dE. All models assume interstellar absorption us¬ 
ing the cross sections of IBalucihska-Church &: McCammonI 1119921b assuming solar 
abundances with h eld fixed at 8 .6 x 10^^ cm~^. We have implemented the 
fast pile-up algorithm of|D avi illool. 


servatories that might arise from the different responses 
of ROSAT and Chandra, the Chandra images are pro¬ 
duced from events restricted to the 0.1-2.5 keV energy 
band of ROSAT. The resultant Chandra images, as well 
as the ROSAT HRI image, are then convolved with a 
Gaussian of FWHM 25" resolution to match the angular 
resolution of the ROSAT PSPC. Finally, to ensure each 
data set has approximately the same number of total 
counts as the ROSAT PSPC and HRI observations, only 
a fraction of the Chandra observations were used to gen¬ 
erate the images as discussed in The images clearly 
show variability in the jet structure and brightness over 
a twelve-year period, supporting the flux measurements 
from the spectral fits. Indeed, the jet is dim and poorly 
defined in 1991/1992. In 199^ however, part of the jet 
structure ~ 3' (corresponding to a distance of 9 pc, cor¬ 
recting for a jet inclination of 30° to the line of sight) 
to the southeast of the pulsar has brightened over the 
two years since the previous observation. By 2000, the 
jet is well-defined with structure extending 4' from the 
pulsar and it has a curved appearance. The 2003 Chan¬ 
dra observation shows that part of the jet ~ V to the 
southeast of the pulsar has dimmed and the structure at 
the end of the jet has become less well-defined. The per¬ 
sistence of the inner- and outer-arc structures during the 
entire time range, taken with variability of the jet struc¬ 
ture and relative intensity, indicates that the observed 
changes are in fact significant and real since it is difficult 
to imagine an artifact that would only apply to discrete 
parts of the images. Hence, the jet brightening indicates 
the existence of an underlying mechanism, as discussed 
in U 

Despite remaining steady in total flux between 2000 
and 2003, the jet shows variability on ^ 20" size scales, 
as shown in Figure^l In 2000, we identify four large (20", 
physical size of 0.5 pc) clumps in the jet (labeled as 1—4 
on the top panel of Figure HI) located between ~ 1' (a 
separation from the pulsar of 3 pc at a 30° inclination) 
and ~ 2.'5 away from the pulsar. In 2003, there are two 


~ 20"-sized clumps located ^ 1.'5 and 1.'7 away from the 
pulsar with fainter emission further down the jet. If we 
interpret the clumps as the same structure having moved 
along the jet, then the velocity is ~ 0.5c (assuming a 30° 
inclination to the line of sight and correcting for rela¬ 
tivistic Doppler boosting). Due to the loss of dithering 
for the 2003 Apr 28 observation, we cannot determine if 
the jet structure changed on week-long time scales. How¬ 
ever, the 10 ks of data in the 2003 Apr 21 observation 
shows jet clumps with similar sizes and locations as in 
the 2003 Oct observation indicating that the jet structure 
changes on timescales longer than 6 months. 

4.5. The Inner and Outer Arcs 

ROSAT images indicated th e presence of a cross-lik e 
structure around the pulsar ijBrazier fc Becked IT9ff7|l. 
With the improved resolution of Chandra, the “cross” 
was res olved into inner and outer arcs (features “5” and 
“E” of IG02I1 located approximately 11” and 40" (dis¬ 
tances of 0.4pc and 1.2 pc at a 60° inclination perpen¬ 
dicular to the jet) away from the pulsar. To determine if 
the arcs have changed structure, we constructed Chan¬ 
dra and ROSAT HRI images matched in energy range, 
counts, and resolutio n to the ROSAT PSPC images as 
described in m and EH Due to the ^ 25" resolution 
of the ROSAT PSPC detector, emission from the pul¬ 
sar and small-scale knots surrounding the pulsar extends 
into the inner and outer arc structures. Therefore, we 
could not exclude the pulsar before smoothing without 
adversely affecting the arc structures. Figure|31shows the 
comparison between the ROSAT and Chandra images. 
There appears to be structural changes to both the inner 
and outer arcs between 1991/1992, 1994, 2000, and 2003. 
In the same manner as for the diffuse PWN in EEI we 
performed simultaneous absorbed power law spectral fits 
to the outer arc using the ROSAT PSPC and Chandra 
data holding Nh fixed at 8.6 x 10^^ cm“^, requiring F to 
be the same for all the data sets, and allowing the normal¬ 
ization to vary between epochs. For the Chandra data. 
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Fig. 3.— Chandra and ROSAT images showing changes in the outer arc and jet from 1991/1992 to 2003. To facilitate compar ison 
between different observatories, the images are matched in resolution (25'^'), energy (0.5—10.0 keV), and counts as described in Oand ^4.41 
The grayscale is linear and ranges from 2.2% to 13% of the peak value for the ROSAT PSPC and Chandra images and 3.9% to 13% of the 
peak value on the ROSAT HRI image to account for the higher background of the HRI detector. 


we reproduce the flux and photo n inde x reported bv lG02l 
within the errors. As noted by IG02L the outer arc has 
the same photon index as the jet and the photon index is 
flatter than for the diffuse PWN. The outer arc appears 
to have decreased in brightness between 1991/1992 and 
2000 by ~20% as shown in Tabled however this is only 
a ItT result (68% confidence limit). This slight flux de¬ 
crease may just be an artifact of the ~ 25" resolution of 
the ROSAT PSPC detector and the pile-up of the pulsar 
in the Chandra data (the low resolution of ROSAT re¬ 
sults in pulsar counts contaminating the outer arc region 
and the pile-up of the pulsar with Chandra results in less 
contamination from the pulsar at the same resolution). 

If we just consider the Chandra data, the outer arc 
shows time variability between 2000 Aug and 2003 Oct 
as shown in Figure 0 In 2000 Aug, the outer arc con¬ 


tained two major clumps - the smaller clump having a 
10" diameter, and the larger clump elongated along the 
arc with a size of 10" x 20" (0.25 pc x 0.5 pc). In 
2003 Oct, the clump locations had changed within the 
outer arc, although the clump sizes were comparable to 
the 2000 clumps. Similar clump variability is also seen in 
the outer arc between 2003 Apr 21 and 2003 Apr 28 as 
shown in Figure IHl indicating that the timescale for arc 
variability might be as short as a few days. The white 
arrows on Figures|31and|niidentify locations of transverse 
structural change between 2000 and 2003. Note that the 
changes observed in the outer arc between the 2003 Apr 
images are only slightly greater than one would expect 
from Poisson statistics. The greatest variation between 
the two 2003 Apr images, indicated by the rightmost ar¬ 
row on FigureEl is only 3cr on unsmoothed images. Spec- 
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Fig. 4.— Chandra images convolved with a Gaussian of FWHM 
10” showing changes in the jet. The numbers indicate the clumps 
as described in the text. The energy range is 0.3—S.OkeV. The 
grayscale is linear and ranges from 0.7% to 7% of the peak value. 


tral fits indicate that the total brightness of the outer arc 
has remained the same from 2000 to 2003. 

The inner arc also shows a brightening to the west of 
the pulsar between 2000 and 2003, as seen in Figures 
01 and 13 however this is most likely due to the bright, 
pointlike source near (J2000) R.A. 15^13™5F;5, deck 
—59°08'25" in the 2003 Oct image seen most clearly in 
the bottom panel of Figure 01 A SIMBAD search shows 
no variable source in this location. While it is possible 
that this source could be related to the small-scale knots 
identified in sa we believe that it is more likely an ar¬ 
tifact of a “hot pixel” in the data. The southwestern 
tip of the inner arc does appear to vary slightly between 
2003 Apr 21 and 2003 Apr 28 (see Figure EJ, however, 
this may be due in part to the lower signal-to-noise ratio 
on the 10 ks exposures. 

The inner and outer arcs were predicte d to h ave out¬ 
ward motions of a few arcseconds per year ll(40l . To de¬ 
termine the proper motion of the outer arc, we convolved 
the 2000 Aug and 2003 Oct Chandra images with a 
Gaussian of FWHM 10". We then constructed an angle- 
averaged radial profile over the entire length of the outer 
arc and centered on the pulsar for each epoch, as shown 
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Fig. 5.— Chandra images convolved with a Gaussian of FWHM 
10'' showing changes in the inner and outer arcs over a three-year 
timescale. The white arrows indicate locations of structural change 
between 2000 and 2003. The energy range is 0.3—S.OkeV and the 
exposure time in each image is 20 ks. The grayscale is linear and 
ranges from 1.4% to 5.9% of the peak value. 
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Fig. 6.— Chandra images convolved with a Gaussian of FWHM 
10" showing changes in the inner and outer arcs over the span 
of a week. The white arrows are the same as in Figure ^ The 
energy range is 0.3—S.OkeV and the exposure time in each image 
is 10 ks. The grayscale is linear and ranges from 1.4% to 5.9% of 
the peak value. The images have been blanked to omit chip gap 
and dead column areas that could not be exposure corrected in the 
2003 Apr 2S observation. 
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Fig. 7. — Angle-averaged radial profiles for the outer arc for 
the 2000 Aug and 2003 Oct Chandra images convolved with a 
Gaussian of FWHM 10''. The difference in brightness at a radius 
of 31" — 32" is due to the variation of small-scale, compact knots 
near the inner arc between 2000 and 2003. Although the small- 
scale knots are unresolved at 10" resolution, they do contribute 
significantly to the fiux of the inner arc region. Possible inward 
motion is indicated. 


in Figure 13 We measured the motion of the brightness 
profiles by minimizing of the difference between the 
profiles at each epoch as a function of radial shift and 
amplitude scaling factor. We used data between radii 
of 38" and 55" (corresponding to the peak and outside 
edge of the brightness profiles) to measure the motion. 
It was necessary to use the peaks of the brightness pro¬ 
files to break the degeneracy between motion and scaling 
differences between epochs. Data at radii less than 38" 
were excluded for the motion measurement because of 
the brightness variations caused by the small-scale knots 
near the inner arc. Although the small-scale knots are 
unresolved at 10" resolution, they do contribute signif¬ 
icantly to the flux of the inner arc region. The error 
was determined at the 68% confidence limit of the 
distribution. We confirmed the error by repeating the 
measurement for ten Monte Carlo simulations of the 
Poisson noise and calculating the rms scatter of the min¬ 
imum x^ position. The outer arc appears to have moved 
inward by 1"0 ± 0'.'2 over the 3-year timespan from 2000 
to 2003. This corresponds to a velocity of 0.03c (assum¬ 
ing outflow perpendicular to the jet with a 60° inclination 
to the line of sight and correcting for relativistic Doppler 
boosting). The convolution of the images and the angle¬ 
averaging of the outer arc results in a measurement of 
the motion of the ensemble average of the material com¬ 
prising the a rc wh ich allows us to make a comparison to 
the model of lC02l Thus, we can report that the average 
arc structure appears to move inward. However, given 
the degree of structural change in this timespan, we can¬ 
not definitively determine if this represents a real inward 
motion or is just a result of clump variation. We can rule 
out unaliased outward motion of the outer arc at the 5cr 
level. For the inner arc, there was too much structural 
change due to the nearby small-scale, compact knots to 
determine proper motions for that feature. 

5. SPATIALLY RESOLVED SPECTROSCOPY 


Although the jet structure has changed over the 3-year 
timespan of the Chandra observations, individual spec¬ 
tral fits to the 2000 Aug and 2003 Oct data indicate that 
there was no photon index evolution in the jet and dif¬ 
fuse PWN during that time as shown in Table Q] To 
determine if the photon index of the diffuse PWN and 
jet varied spatially, we used just the four Chandra data 
sets and performed simultaneous absorbed power law fits 
as outlined in ED For each region described below, we 
held Nh fixed to 8.6x 10^^ cm“^, required the same F for 
all the data sets, but allowed the normalization to vary 
between epochs. 

For the diffuse PWN, we extracted spectra from con¬ 
centric 36"-wide annuli extending from a radius of 113 
to 4.'3, centered on the pulsar as outlined in Table 0 
The background region used for the diffuse PWN was a 
36" wide annulus just exterior to the outermost “source” 
annulus and interior to the RCW 89 region. The back¬ 
ground annulus was exterior to most of the PWN ex¬ 
cept to the northwest and to the southeast where it in¬ 
cludes faint PWN emission. The jet was excluded from 
the diffuse PWN analysis. We performed a linear least- 
squares fit to determine the rate at which the photon 
index changed with radius. The results are shown graphi¬ 
cally in FigurelHl The diffuse PWN shows a steepening of 
0.04±0.02 arcmin“^ between 113 and 3ll (for a total AF 
of ~ 0.07) but flattens significantly thereafter. If we al¬ 
low Nh to also vary for the diffuse PWN fits, the photon 
index does not flatten at large radii and the best fit slope 
is 0.08±0.02 arcmin”^ which translates to a total change 
in photon index between ll3 and 414 of 0.25. Nh in¬ 
creases from 9.2x10^^ cm“^ in the annulus closest to the 
pulsar to 13.2x10^^ cm“^ in the annulus furthest from 
the pulsar. These values for the absorbing column are 
consistent with those derived by iTrussoni et al ] iTT^ 
showing that there was about 1.5 times more absorption 
towards RCW 89 than towards the pulsar. 

For the jet, we constructed five regions along the length 
of the jet as shown in Figure ^ The separation from the 
pulsar was defined at the center of each region. The 
background for the jet analysis was chosen as the region 
immediately exterior to the total jet region. A linear 
least-squares fit was applied to the data and is shown 
in Figure |3 The jet shows a marginal photon index 
steepening of 0.17±0.09 arcmin“^ between radii of 1' to 
315 for a total change of 0.43. 

6. DISCUSSION 
6.1. Variability of the Arcs 

IG02I suggested that the inner and outer arcs represent 
compressions downstream of the termination shock in an 
equatorial electron-positron pair outflow. The compres¬ 
sions are induced as the high energy heavy ions (“pro¬ 
tons”) embedded in the flow enter the shock heated 
downstream pairs with the same velocity and Lorentz 
factor as the upstream pair wi nd. This was an applica- 
tion of the model advanced bv iCallant fc Aron^ lll994ll 
for the time variable wisps seen in the termination shock 
region of the Crab Nebula , which has been shown by 
iSpitkovskv fc Aroii^ ll2004 hereafter SA04) to give a 
good representation of the structure and toroidally av¬ 
eraged variations seen in the inner X-ray ring as well as 
the propagating wave structures eme rging from that ring 
seen by Chandra llHester et al.ll20d^ . In the application 
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TABLE 3 

Radial variations in photon index for the 

DIFFUSE PWN AND JET. 


Region 

Distance 

(') 

Njf 

(10^^ cm“'^) 

F 

Jet'*' 

1.0 

8.6 (fixed) 

1.48±0.27 


1.5 


1.49±0.26 


2.0 


1.70±0.16 


2.7 


1.86±0.24 


3.5 


1.64±0.45 

Diffuse EWN*" 

1.3 

8.6 (fixed) 

1.84±0.03 


1.9 


1.91±0.03 


2.5 


1.88±0.03 


3.1 


1.93±0.02 


3.7 


1.84±0.03 


4.3 


1.74±0.08 

Diffuse PWN 

1.3 

9.2±0.5 

1.91±0.05 


1.9 

9.7±0.5 

2.02±0.05 


2.5 

10.3±0.5 

2.05±0.05 


3.1 

10.6±0.6 

2.14±0.07 


3.7 

10.9±0.8 

2.07±0.09 


4.3 

13.2±1.4 

2.18±0.14 


^The regions used to extract spectra for the jet are 
shown in Figure 13 Distances from the pulsar are mea¬ 
sured at the geometric center of each region. 

^The regions used to extract spectra for the diffuse 
PWN are concentric 36^^ wide annuli. Radius is given 
at the center of each annulus. The jet region has been 
excluded from the diffuse PWN analysis. 


of that model to G320.4—1.2, the inner and outer arcs 
were interpreted as the turning points in the ions’ orbits 
as they do their first gyrations in the sho ck com pressed 
magnetic field. Fitting the model to the lGfl2l observa¬ 
tions of G320.4—1.2 led to an estimate of the pair flux, 
upstream flow Lorentz factor, ion flux a nd magnetic field 
in the equatorial flow. From that. iGf)^ derived an ion cy¬ 
clotron time on the order of one to two years, which is the 
basic variability time of the ion induced compressions. If 
all of the variation is in the f orm o f bulk motion, the 
predicted outflow velocity from iGf)^ is ^ 0.5c. Although 
our observations do not show an outward velocity, we 
do confirm the theoretical prediction of the variability 
time scale expected in the inner and outer arcs. That 
model, being a toroidal average, did not predict the an¬ 
gular structure of the arcs uncovered by the new data 
presented here. 

This model depends on the existence of ultra high en¬ 
ergy heavy ions being present in the equatorial flow, with 
the ions carrying a large enough fraction of the energy 
flux to be able to induce substantial compression in the 
radiating pair plasma. The inferred ion number flux, 
both in the Crab and in G320.4—1.2, is approximately 
what one expects on electrodynamic grounds, if the ion 
stream is the electric return current required to prevent 
the charging up of the neutron stars. It is not cur¬ 
rently understood how the particles in the flow achieve 
their high energies (7 ~ 10 ®), although it has been sug¬ 
gested that they are accelerated and heated by the mys¬ 
terious dissipation processes that lead to the equatorial 
wind h aving magnetization as weak as has been inferred 
lISAOl and references therein). 

These kinetic models focus on the time dependent 
morphology in and around the termination shock in 
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PWN with linear least-squares fits (dashed lines) and the slope 
of the fit indicated. Error bars are shown for the 90% confidence 
limits. 
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the pairs - in effect, the ions form part of a re¬ 
solved shock structure whose basic form is an equa¬ 
torial ring. The large scale structure of the nebular 
flow has been addressed in a series of MHD models 
and simulations, whose focus is on understanding the 
dramatic pole to equa tor asymmetry, especially the ap¬ 


pearance of the jets llBogovalov fc Khang qulianl 
Lviibarsl^ 120021: IKomissarov fc Lviibarsk-yl 120031 


2002 ; 


2004; 


De^anna^^nat(^7^Bucciantmil20(3Mogo^!avm^Lr] 

2005^ . These models assume all the particles have 
small Larmor radius, and thus are indifferent to com- 
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position and exclude the proposed high energy ions by 
assumption. With the premise that the injected en¬ 
ergy flux decreases with increasing rotational latitude 
Q (usually assumed proportional to cos^ Q ), the energy 
is mostly injected through a toroidal belt shock, with 
some flow passing obliquely through a higher latitude 
“arch shock” which, because the shock is oblique to the 
flow, does not decelerate the plasma to subsonic speeds. 
These non-spherically symmetric shock structures can 
create the appearance of an inner ring of emission, with 
the weakly decelerated plasma emerging from the arch 
shocks contributing to form the appearance of an outer 
ring. Doppler boosting may yield an appearance of par¬ 
tial arcs rather than complete rings®. These models pro¬ 
vide a plausible scenario for the formation of the jets, as 
the consequence of backflow and magnetic hoop stress at 
higher latitudes. 

The simulations of this model show no signs of short 
time variability corresponding to the wisps®. The simu¬ 
lations may lack suflicient resolution to And the Kelvin- 
Helmholtz instability between the outward flow emerging 
from the arch shock and the hi g her la titude backflow, an 
effect suggested by iBegelmaiil lll999|l with a somewhat 
different flow geometry in mind as the origin of the wisps. 
Alternatively, magnetization high enough to form the jet 
may suppress the shear flow instabilities. 

While the effects of high energy equatorial ions have 
not been explicitly included in the MHD simulations, if 
present they would be injected into the shock decelerated 
pairs through the equatorial ring shock seen in the MHD 
simulations. Thus the formation of ion driven compres¬ 
sions is qualitatively consistent with the MHD models, 
and remains the only quantitatively elaborated model for 
the time variability seen in the Crab wisps and perhaps 
seen in the data described here. Based on this model, 
IC02I predicted arc motions of a few arcseconds per year. 
We see no outward moti on of t he outer arc and perhaps 
even an inward motion. ISA04I noted that the positions 
of the ion driven compressions can move inward as well 
as outwards, with the apparent motion then depending 
on when one takes a snapshot of the structure - the 
possible existence of apparent inward motion, perhaps 
seen in our data, is a prediction of the model. However, 
the apparent motion may also be due to the transverse 
structural variations in the outer arc between 2000 and 
2003 - non-axisymmetric variability, a 3D effect, has not 
been included in any of the models published to date. 
If such transverse change is the origin of the apparent 
motion, the arcs would then be quasi-stationary in the 
ra dial sens e . Or, we could be witnessing alaising as seen 
bv lSca.rgl^ ljl969jl in his observations of the wisps in the 

® However, such a model then has difficulty explaining why the 
inner X-ray ring in the Crab Nebula has an approximately axisy- 
metric appearance, after toroidally averaging over the knots in the 
ring brightness, while the larger radius torus does show signs of 
Doppler boost. The ion compression model has the same difficulty. 
One of the major successes of the MHD model is to find a flow that 
does have substantial Doppler boost at distances from the pulsar 
comparable to the Crab torus 

® |Bogovalov_et.^f3 (2Q23) do report a long period quasi-coherent 
oscillation of the flow (t ~ 40 years) with some possible shorter 
time variability when the magnetization is finite but not large - the 
variability vanishes for unmagnetized flow models. The physical 
reasons for this variability are not understood (Bogovalov, private 
communication). 


TABLE 4 

Alfven Crossing Times in the Jet 


Position^ 

-®min 

ifiG) 

^min 

erg cm~^) 

(years) 

1 

15 

2.1 

1.3 

2 

17 

2.8 

1.3 

3 

14 

1.7 

2.0 

4 

12 

1.3 

2.0 

5 

11 

1.1 

1.6 


'‘Positions are shown in Figure Q position 1 is 
closest to the pulsar. Volumes at each position are 
assumed to be cylindrical. 

'‘The Alfven velocity at each position along the 
jet as computed from Equation^is ^ 0.63c. 

Crab, due to the observations having undersampled the 
variations. 

We can test whether the quasi-stationary nature of the 
outer arc is due to aliasing. If we assume that the par¬ 
ticles comprising the outer arc are relativistic, and thus 
Doppler boosted, then we can determine a speed based 
on the brightness ratio of > 5 between the near and far 
sides. For a photon index of 1.6, the Doppler boosting 
formula reduces to Pcoscj) = 0.22, where jS = v/c and (f) 
is the inclination angle to the line-of-sight. If the incli¬ 
nation of the outer arc is 60° (oriented 90° from the jet), 
then (3 = 0.44. At this speed, we expect the outer arc to 
move nearly 16" in 3 years, which is slightly greater than 
the width of the arc. Over 6 months, the outer arc should 
move 2" — 3", which is not observed. Suppose that the 
motion of particles through the arc is not indicative of 
the motion of the arc. In this case consider that the arc 
would need to move at 0.27c to move its width (0.25 pc) 
in 3 years. At this speed, the arc would only move 176 in 
6 months. This amount of motion is difhcult to rule in 
or out given the degree of structural change in the outer 
arc and the signal-to-noise ratio of the data. 

The nature of arc a and arc b is unknown. If indeed the 
arcs are outward-moving structures, then perhaps arc a 
and arc b are just more evolved versions of the inner and 
outer arcs. They could also represent a quasi-coherent os¬ 
cillation of th e flow in the outer PWN as observed in the 
simulations of iBogovalov et ahl (120051) . A more thorough 
analysis of these faint PWN structures will be performed 
at a later time. 

6 .2. Variability in the Jet 

MHD sausage or kink instabilities are via ble candidates 
for th e jet variability, as in the Vela jet llPavlov et alJ 
Eqq 3). The time scale for such phenomena is set by the 
Alfven crossing time tA = rjet /f a, where r^et is the cylin¬ 
drical radius of the jet and 

= , « , ^ ( 1 ) 

Here p is the relativistic pressure, assumed to be 
isotropic, and U = 3p ^ pc^ is the energy density of 
the relativistic plasma. Shown in Table 0 are the mag¬ 
netic fields, plasma energy densities, and Alfven crossing 
times for regions along the jet as shown in Figure Q] Po¬ 
sition 1 is closest to the pulsar and we are assuming a 
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cylindrical volume. The magnetic fields and energy den¬ 
sities are computed from equipartition arguments using 
integrated fluxes for the energy range 0.5-10 keV and 
assuming an uncooled spectrum llPacholczvld 1197(1 pp. 
170-171). For all regions va ~ 0.63c and the Alfven 
crossing times vary from 1.3—2 years. Thus the MHD 
instability of a magnetic pinch, the basic structure of the 
jet described in the MHD models, certainly is a viable 
candidate for the jet variations observed between 2000 
and 2003, and perhaps may be a candidate to explain 
the partial brightening of the jet from 1991/1992 to 1994 
and 2000. Indeed, episodic jet outflows are observed in 
MHD simulations for certain magne t ic field configura- 
tions llOuved. Pudritz fc Stond Il997t lOuved fc Pudrit j 
[1993). 

If we interpret the jet structural changes between 2000 
and 2003 as clump motion along the jet, then the result¬ 
ing velocity is ~ 0.5c. This velocity is consistent with 
the prediction of IG02I based on Doppler boosting and 
also matches the outflow speeds obs erved in the Crab 
Nebula jet (0.4Cj_|Hester_e£MTj20^, the Vela PWN 
jet (0.3c — 0.7c, iPavloy e t alJ 1200^)7 and in the jet of 
Gil.2—0.3 (0.8c—1.4c. iRoberts et al.ll2003j) . Relativistic 
jet outflow velocities of ~ 0.5c are predicted by the 
MHD si mulations described above llBogovalov et al.l 

2^2^ _ iDel Zanna. Arnato. fc Buc . ciantinil 12004 

Komissarov fc Lvubarskv I2003L l2004|) . The jet bright- 
nening from 1991/1992 to 1994 and 2000 cannot simply 
be a result of material moving at 0.5c down the jet 
causing the jet to lengthen and brighten with the 
inclusion of more material. This scenario is excluded 
for two reasons, the first is that it would take nearly 80 
years for material moving at 0.5c to travel the length 
of the jet, and the second is that the jet seems to 
brighten in a “patchy” manner and does not simply 
lengthen, as shown in Figure O Therefore, as discussed 
above, the brightening of the jet is more likely the 
result of an MHD instability such as a magnetic pinch. 
In this model, the instability causes compression and 
brightening of local regions, each of which has length 
no more than the jet diameter. There is no necessity 
for the energy to be injected into the jet at a varying 
rate, with the observed variations reflecting transport 
along the jet. This scenario works because the Alfven 
transit time across the jet is much shorter than the flow 
time along the jet. The model suggests that the clumps 
vary incoherently with respect to each other, on a time 
scale approximately equal to the Alfven crossing time 
(< 2 years), consistent with the very limited time series 
available to us. 

6.3. Small-Scale Structure Near the Pulsar 

The compact, small-scale knots near the pulsar are 
quite variable on 6 month timescales. Although the knot 
to the southeast of the pulsar might represent material 
moving at 0.6c, given the startling knot variability to 
the northwest, a conclusion of knot motion is perhaps 
premature. We can estimate the magnetic fields in the 
knots in three ways: 1. assume that the knot size is ap¬ 
proximately the Larmor radius of gyrating pairs, 2. use 
equipartition/minimum energy, and 3. use synchrotron 
lifetime if the knots are in motion at 0.6c and cool be¬ 
low detectability at a distance of the inner arc in time 
frames from 6 months to 1 year. For the first method, 


IG02I derived a magnetic field of 3/rG from the Larmor 
radius for their knot 1. For the second method, the mag¬ 
netic field derived from minimum energy us ing a size of 
0.1 pc and the integrated flux of knot 1 from IG03 for the 
energy range of 0.5—10.0 keV is Hmin = 76^G. For the 
third method, the synchr otron lifetime can be calculated 
using (from equation 4 of lG02ll 

tsynch = kyr (2) 

where B is the magnetic field in faG and e is the energy 
in keV. For 5 keV and tsynch = 1 year, B = 700/iG. 
Given this high magnetic field, it seems unlikely that 
synchrotron cooling could be responsible for the transient 
nature of the knots. 

The knots may have a similar explanation as “knot 2” 
in the Crab Nebula which has been interpreted as an 
unstable quasi-stat ionary shock in the polar jet outflow 
llHester et al.ll20(i^ . In this scenario, one would expect 
the bulk of the knots to be in the approaching jet since 
the receding jet is Doppler boosted down in brightness. 
However, in this case we observe most of the knot struc¬ 
ture on the receding jet side. 

Given that the knots are never observed beyond the in¬ 
ner arc radius and most of them occur on the approach¬ 
ing side of the arcs, it is likely that the knots are ei¬ 
ther associated with the equatorial outflow or the result- 
in g backflow as described in il6TJand shown in Figure 3 
of iKomissarov fc Lvubarsl^ ll2003tl . Two processes ex¬ 
pected from the relativistic MHD simulations that could 
produce “knot-like” emission are Kelvin-Helmholtz (KH) 
instabilities arising in the shear layer between the equato¬ 
rial outflow and the backflow and unsteady vortices that 
develop in the converging flow at the base of the jet. 

The KH instability creates waves that may grow and 
break, while traveling in the direction of the mass 
weighted jump in velocity across the interface between 
the outflow and backflow. Doppler boosting as the break¬ 
ing waves head toward us, and plasma and B field com¬ 
pression as the waves break, may lead to the formation of 
bright features that can appear to us as localized bright 
spots of emission. The main difficulty with this idea is 
that within the flow models, the mass weighted velocity 
jump is directed toward us only on the near (visible jet) 
side except very close in to the base of the jets. This 
is because flow convergence increases the mass density 
in the layer that converges on the axis at small distance 
from the axis, while the outflowing layer has low density 
at small radii. At larger radii, the relative densities of the 
layers reverses. Thus this effect will work only if there is 
a lot of flow convergence. 

The converging flow at the base of the jet could also 
produce transient features due to the development of un¬ 
stable vortices (N. Bucciantini, private communication). 
The backflow would need to focus quite well toward the 
base of the jet to produce Doppler boosting of flows di¬ 
rected towards us and may be rather sensitive to the 
maintenance of axisymmetry. However, this hypothesis 
naturally explains the relative absence of knots on the 
near (visible) jet side since the converging flow at the 
base of the approaching jet would be directed away from 
us. 

Regardless of the actual cause of the knots, they seem 
to be diagnostic of flow details in the PWN. To determine 
the true nature of the knots will require much higher res- 
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olution simulations and more analysis of the flow models. 

6.4. Synchrotron Cooling in the Jet and Diffuse PWN 

IG02I determined that the synchrotron cooling time for 
the jet at 5 keV and a magnetic field of ~ 25 /rG is 
140 years. This cooling should result in a AF of 0.5 

is a marginal indication of cooling in the jet with a to¬ 
tal AF of 0.43. We await analysis of the new, longer 
Chandra observation to verify the possible synchrotron 
cooling in the jet. The relativistic MHD simulations do 
predict that the jet should undergo synchrotron cooling, 
although there is mixing with the ambient PWN material 
(S. Komissarov, private communication). At the base of 
the jet, the photon index should match that of the equa¬ 
torial outflow since the jet is constructed from purloined 
equatorial material. If we assume that the inner and 
outer arcs identify the equatorial outflow material, then 
their photon index of 1.6 does nearly match the photon 
index of 1.5 at the base of the jet. If the spectral steep¬ 
ening along the jet is real, it would be sim ilar to that 
observed in the Grab Nebula jet lAF ^ 0.5; iMori et al.| 
1 ^ . 


Spectral steepening with increasing radius has been ob¬ 
served in a number of PWNe and is generally attributed 
to both synchrotron cooling and expansion losses 


lIGonzalez h Safi-HarbI 120011 

: ILu et alJ 120021 

Mori et alJ 

120041 ISafi-Harb et a,lJl2001l 

Slane et all 12004 

). The de- 


gree of steepening varies from AF ^ 0.3 to 1.3 with a 
typical value of 1 and the radial photon index profiles 
generally show a monotonic increase with radius. Our 
spectral results for the diffuse component of the PWN 
of G320.4—1.2 also show a monotonic increase of pho¬ 
ton index with radius if we allow the absorption to vary, 
but the degree of steepening (~ 0.25) is less than in 
other PWNe. In contrast, m odels for synchrotron cool- 
ing in PWNe based on the iKennel fc Goronitil lll984f) 
model predict a rapid i ncrease in photon index at large 
radius llRevnoldsll2003t ISlane et al.ll2004ll . Perhaps the 
discrepancy between the cooling model and the observa¬ 
tions is due to large-scale mixing of recently accelerated 
material with “older” material as predicted by the rel¬ 
ativistic MHD models llB ogovalov fc Kha ngoulian] 2002; 
L^ubarsk^j200^_|Komissarov&L^ubarekj|j200^ 2004' 

De^aima^Amatcx"^r^ucciannnl2()0^lBogo^dav^t'^d] 

200 ^ 

7. SUMMARY AND CONCLUSIONS 

Variability is observed in the X-ray PWN of PSR 
B1509—58 on time scales possibly as short as one week 
and up to twelve years. Our primary results are as fol¬ 
lows: 

1. The compact, small-scale knots appearing within 
20 " of the pulsar exhibit transient behavior which 
may be attributed to turbulence in the flows sur¬ 
rounding the pulsar. Possible knot motion is indi¬ 
cated with a velocity of 0.6c. 

2. Apparent outflow along the jet is observed with 
velocities of ~ 0.5c. This outflow alone cannot ac¬ 
count for the ~30% brightening of the jet between 
1991/1992 and 2000. The Alfven crossing time 


for the jet is 1.3—2 years, therefore, MHD kink 
or sausage instabilities can account for the rapid 
morphological variations and perhaps the partial 
jet brightening. 

3. The outer arc has possibly moved inward with a 
velocity of 0.03c, however the transverse structural 
changes seen in the outer arc may account for the 
apparent motion. We cannot determine at this 
time if the outer arc is truly quasi-stationary or 
if we are witnessing aliasing. 

4. The diffuse PWN has not evolved significantly in 
structure or brightness over the 12-year time span. 
Using the summed Chandra images, we identify 
two possible arc structures exterior to the outer 
arc. 

5. The photon indices of the diffuse PWN and possi¬ 
bly the jet steepen with increasing radius indicating 
synchrotron cooling at X-ray energies. 

Although our imaging capabilities have improved sub¬ 
stantially since the first optic al observations of time vari¬ 
ability in the Crab Nebula llScarglelll969ll . our under¬ 
standing of these variations in PWNe is still quite lim¬ 
ited. For instance, while we expect magnetic fields to 
play an important dynamical role in jets, and indeed 
we do see variations on the appropriate Alfven crossing 
times, we do not know for certain if MHD instabilities 
are the root cause of the observed variations. The arc 
structures we observe in G320.4—1.2 are equally enig¬ 
matic. We do not yet know if they are in steady motion 
or are quasi-stationary wave phenomena. The striking 
changes in the small-scale knots near the pulsar may sim¬ 
ply be “weather,” diagnosing unimportant details in the 
PWN flow, or they may indicate important flow structure 
which is essential to understanding, for instance, diffu¬ 
sion of particles from the equatorial flow to higher lati¬ 
tudes, a loss essential to a post pair shock second-order 
Fermi acceleration model. Certainly deeper, and appro¬ 
priately spaced, X-ray observations will help resolve some 
issues such as the possible spatial aliasing of the outer 
arc. Also, the higher signal-to-noise will provide better 
constraints on the spatial spectral index variations and 
allow us to determine if and how much mixing has oc¬ 
curred in the diffuse PWN. Finally, we are excited by 
the recent development of relativistic MHD models and 
we hope that some of the variability we observe here can 
eventually be observed in those simulations. 
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